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In the lattice formulation of the Heavy Quark Effective Theory (LHQET), the "classical velocity" v becomes 
renormalized. The origin of this renormalization is the reduction of Lorentz (or 0(4)) invariance to (hyper)cubic 
invariance. The renormalization is finite and depends on the form of the decretization of the reduced heavy quark 
Dirac equation. For the Forward Time - Centered Space discretization, the renormalization is computed both 
perturbatively, to one loop, and non-perturbatively using two ensembles of lattices, one at P = 5.7 and the other 
at P = 6.1 The estimates agree, and indicate that for small classical velocities, v is reduced by about 25-30%. 



1. INTRODUCTION 

1 9 

In the heavy quark limit ' a single form factor, 
the Isgur-Wise universal function ^, describes all 
semileptonic decays of one meson containing a heavy 

quark into another, such as the process B ->D * / v . 
Its calculation is well suited to lattice gauge 
theory , and such calculations have been carried out 
using a lattice implementation of the heavy quark 
effective theory^, also by treating the heavy quarks 
as Wilson fermions with a small hopping constant, 
without implementing the heavy quark limit^'^. 

On the lattice^, as in the continuum^'^'^, the 
Isgur-Wise limit entails the introduction of a 
"classical velocity", v, normalized to 1, which 
appears both in the decomposition of the momentum 
of a heavy particle and the reduced Dirac equation of 
the heavy quark field. 

P = Mv + p 

(1) 

-ivDh^''\x) = 

In the continuum, the second relation is derived from 
the first, and the velocity that appears in these two 
contexts is the same. However, on the lattice this is 
not the case. 

In this talk, we explain that the origin of this 
difference is the reduction of the space-time 
symmetry group from Lorentz (or Euclidean 0(4)) 



invariance to hypercubic lattice symmetry. We then 
describe two calculations of its magnitude. The first 
is a one-loop perturbative calculation, which follows 
the analysis of Aglietti*^ (who however used a 
different discretization of the lattice Dirac operator 
than that adopted here). The second calculation is 
non-perturbative. It is based on a computation of the 
shift in the energy of a meson containing a heavy 
quark, as measured by the fall-off of its Euclidean 
space propagator, for a given shift in its residual 
momentum. In both calculations, it is useful to 
expand the shift in v in a power series in the "bare" 
classical velocity, and in the conclusions we compare 
the first three terms calculated using each procedure. 

2. ORIGIN OF THE RENORMALIZATION 

The heavy quark effective theory is formulated^ by 
factoring the M— ^ singular behavior from the field 
of a heavy quark, leaving a reduced operator h . 

The Lagrangian for h is 

^^'^ - h^^Kx) ivD h^'\x) 

and its propagator satisfies the reduced Dirac 
equation 



Here D is the covariant derivative. This structure is 
the same both on the lattice and in the continuum. 
In the continuum, the propagator of the reduced 

field h^^Hx) is 



1 



VP - ^^"hp) 

where 'L^^\p) is the proper self mass. In 
perturbation theory, it is the sum of all the 1 -particle 
irreducible self-mass diagrams. 

If 'L^^\p) is Taylor expanded in the residual 
momentum. 



(5) 



m + 



p = 



(6) 



the first term gives a mass shift. This is without 
physical significance and can be eliminated by a 
redefinition of the reduced heavy quark field by a 
phase which is independent of the heavy mass^. The 
renormalized classical velocity is inferred from the 
second term in the expansion. For small residual 
momentum, the full quark propagator has the 
behavior 



(ren)' 



= Z^^> (V, 



= - 1 



(10) 



On the lattice, p is replaced by an appropriate 
discretization. This is the starting point of the 
calculation in perturbation theory. 

The only 4-vector on which can depend is the 
bare classical velocity. In the continuum, the only 
possible dependence is for is to be linearly 
proportional to v^^. Therefore, it modifies the 
inverse of the heavy quark propagator to first order 
in p only by an overall multipUcative factor, which 

is a wave function renormalization of h^^\x), but 
not a change in the classical velocity. 

On the lattice, because of the reduced Lorentz or 
rotational symmetry, need not be linearly 

proportional to v^. For example, it can have a term 
3 

proportional to , which has the same hypercubic 
group transformation properties as v^^, but is linearly 
independent of it. When the coefficient X^ is not 
simply proportional to v^^, the classical velocity is 
shifted. 



3. NON-PERTURBATIVE SHIFT 



(p) = 



.(V) 



m + v^'^^"^-p + 



(7) 



This identifies the physical classical velocity as 

..(ren) ^ Z^'Hv- XJ («) 



where 



A procedure for calculating the renormalized 
classical velocity which is suited to non-perturbative 
simulation is based on the observation that the rate of 
fall-off of the reduced propagator in Euclidean time 
for fixed residual spacial momentum is given by the 
residual energy at that momentum. 



Sit,p) ~ Z^^'e 



(11) 



32: 



(9) 



The wave function renormalization constant Z*-^'^ is 
determined by the requirement that both the bare and 
the renormaUzed classical velocities have unit 
normalization. 



The residual energy, which is directly amenable to 
simulation, has a simple dependence on the physical 
classical velocity and the residual momentum. In the 
M—^ limit, the shift in the energy is 



E^'\p) 



(12) 



where the conventional classical energy (which is 
bounded by 1) is 




quark proper self mass one loop proper self mass 



; (phys) 
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(13) 



For both the calculations, we employ the 
forward-time symmetrical- space discretization of the 
reduced Dirac equation 

[U{x,x+t)S{xH,y)-S{x,y)] 
3 

+ E ^ [Uix,x^(i) Six^(i,y) ^^4^ 
-U{x,x-{i) S{x-ii,y) ] 

= 8(x,y) 

The use of an asymmetric time difference is required 
in order to implement the requirement that heavy 
quarks propagate only forward in time. 

The free inverse propagator corresponding to this 
discretization is 



,(0)-l 



ip ^ (15) 

M = l 



As is evident from Eq. (14), the simulated 
propagator is naturally expanded in a power series in 
the (bare) classical velocity, which plays the role of 
a transverse hopping constant. To facilitate 
comparison of the perturbative and lattice 
calculations, we will also express the results of the 
perturbative calculation in powers of the bare 
classical velocity. 

3. PERTURBATIVE EVALUATION 

To one loop, the proper self mass is given by the 
two diagrams shown in Figure 1. The point 



interaction only gives rise to a residual mass, and so 
can be ignored. The second diagram does change 
the classical velocity, however. Because of the finite 
lattice spacing, the integration domain is periodic in 
the Euclidean 4-momentum. The heavy quark 
propagator is as in Eq. (15). 

The evaluation of the loop integral over the gluon 
4-momentum requires that the contour in the 
complex Euclidean energy plane be chosen so that 
the heavy quark propagator always vanishes for 
negative Euclidean time. This requires that in the z 
plane it encloses the pole coming from the heavy 
quark propagator, as is indicated in Figure 2. 

The lowest order perturbative renormalization 



o - - (ren) 

dV; = V; - V; 



has the structure (C2 is the quadratic Casimir) 
The expansion coefficients are evaluated numerically: 



'12 



.17800996 
.03159228 
.02959309 



(18) 



4. SIMULATION 

To avoid the ambiguities associated with global 
gauge fixing, we simulate the propagator of a heavy- 
light meson. Its classical velocity is the same as its 
heavy quark component. We have simulated the 
propagator expansion coefficients on an ensemble of 



lattices and Wilson light quark propagators made 
available to us by the Fermilab ACP-MAPS 
Collaboration'*^. These consisted of 48 lattices of 
size 24^ X 48 with lattice coupling P = 6.1 along 
with Wilson quark propagators with hopping constant 
K= .154. We computed the coefficients of the 
leading terms in the expansion of the heavy quark 
propagator, and formed heavy-light meson 
propagators expansion coefficients. 

In this preliminary analysis, we did not optimize 
the wave function of the heavy-hght meson, and in 
fact there was no extended region in Euclidean time 
over which the ground state was cleanly isolated. 
With that caveat, the propagator terms from which 
the physical classical velocity was extracted is shown 
in Figure 3. 

5. COMPAMSON OF RESULTS 
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Figure 3 — Coefficients of the leading terms in 
the logarithmic derivative of S*^^^ 



The perturbative and simulation calculations of 
the coefficients of the first three terms in the bare 
classical velocity expansion of the finite 
renormalization of classical velocity is shown in the 
following Table: 



Table I — Comparison of Simulated and 
Perturbative Expansion Coefficients 



Term 


Simulation 


Naive 
Perturb 


Tadpole 
Improv 




-.287 + .092 


-.2334 


-.3969 




-.982 + .087 


-.0414 


-.0704 




-.764 + .078 


-.0461 


-.0783 



The coefficient of the linear term agrees rather 
well between the simulated and perturbative 
evaluations. There is a statistically insignificant 
preference for the use of the naive rather than the 
tadpole improved coupling in this context. However, 
the perturbative and simulated cubic coefficients are 
in complete disagreement While it is possible that 
the lattice coupling is simply too large to neglect 
higher order terms, is should be pointed out that each 
of the cubic coefficients is the difference of two 
terms, each of which is substantially larger than their 
difference. Furthermore, each of the two terms has 
a t ^ asymptotic behavior, while their difference 



should only grow linearly. The large cancellations 
are roughly of the same magnitude in both cubic 
coefficients. Fortunately, the dominant, linear 
coefficient is free of this difficulty. 

Both the simulation and perturbation theory 
indicate that the physical classical velocity is smaller 
than its bare value. 

REFERENCES 

1. M.B. Voloshin and M.A. Shifman, Yad. Fiz. 45, 
463 (1987); Sov. J. Nucl. Phys. 47, 511 (1988). 

2. N. Isgur and M.B. Wise, Phys. Lett. B232, 113 
(1989); B237, 527 (1990). 

3. C/M. Neubert, Physics Reports 245, 259 (1994). 

4. J.E. Mandula and M.C. Ogilvie, Phys. Rev. D 
(Rap. Comm.) 45, R2183 (1992); Nucl. Phys. B 
(Proc. Suppl.) 26, 459 (1992); 34, 480 (1994). 

5. C. Bernard, Y. Shen, and A. Soni, Phys. Lett. B 
317, 164 (1993); Nucl. Phys. B (Proc. Suppl.) 
30, 473 (1993); 34, 483 (1994). 

6. S. P. Booth et. al. (UKQCD Collaboration), 
Phys. Rev. Lett. 72, 462 (1994); Nucl. Phys. B 
(Proc. Suppl.) 34, 471 (1994). 

7. H. Georgi, Phys. Lett. B 240, 447 (1990). 

8. U. Aglietti, HEP-PH/9304274 (1993). 

9. A. F. Falk, M. Neubert, and M. E. Luke, Nucl. 
Physics B 388, 363 (1992). 

10. See, for example, E. Eichten et al., Nucl. Phys. 
B (Proc. Suppl.) 34, 444 (1994). 



